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Abstract 


Vero E6 African green monkey kidney cells are highly susceptible to infection with the newly emerging severe acute respiratory syndrome 
coronavirus (SARS-CoV), and they are permissive for rapid viral replication, with resultant cytopathic effects. We employed cDNA microar- 
ray analysis to characterize the cellular transcriptional responses of homologous human genes at 12 h post-infection. Seventy mRNA tran- 
scripts belonging to various functional classes exhibited significant alterations in gene expression. There was considerable induction of heat 
shock proteins that are crucial to the immune response mechanism. Modified levels of several transcripts involved in pro-inflammatory and 
anti-inflammatory processes exemplified the balance between opposing forces during SARS pathogenesis. Other genes displaying altered 
transcription included those associated with host translation, cellular metabolism, cell cycle, signal transduction, transcriptional regulation, 
protein trafficking, protein modulators, and cytoskeletal proteins. Alterations in the levels of several novel transcripts encoding hypothetical 
proteins and expressed sequence tags were also identified. In addition, transcription of apoptosis-related genes DENN and hIAP1 was upregu- 
lated in contrast to FAIM. Elevated Mx1 expression signified a strong host response to mediate antiviral resistance. Also expressed in infected 
cells was the C-terminal alternative splice variant of the p53 tumor suppressor gene encoding a modified truncated protein that can influence 
the activity of wild-type p53. We observed the interplay between various mechanisms to favor virus multiplication before full-blown apoptosis 
and the triggering of several pathways in host cells in an attempt to eliminate the pathogen. Microarray analysis identifies the critical host— 
pathogen interactions during SARS-CoV infection and provides new insights into the pathophysiology of SARS. 
© 2005 Elsevier SAS. All rights reserved. 
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1. Introduction 


In November 2002, a novel infectious agent causing severe 
acute respiratory syndrome (SARS) emerged in Guangdong 
province, China. However, a global outbreak alert was 
declared by the World Health Organization (WHO) only in 
March 2003. By July 2003, SARS had spread to about 
30 nations, resulting in more than 8000 cases with either atypi- 
cal pneumonia or respiratory distress syndrome. When WHO 
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finally declared the end of the global outbreak in July 2003, 
SARS had already claimed more than 770 lives. In Singapore 
alone, there were 238 cases with 205 patients who survived, 
but 33 fatalities, accounting for a mortality rate of approxi- 
mately 14%. 

During March 2003, a number of laboratories indepen- 
dently reported the isolation of a novel coronavirus (CoV) 
from SARS patients [1-3]. Thin-section electron micros- 
copy revealed viral particles similar to the morphology of 
CoV. Sequence analyses demonstrated significant similarity 
of SARS-CoV to the family Coronaviridae, especially with 
members belonging to group II CoVs, although the novel CoV 
was genetically distinct from other known CoVs. 
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There have been reports on clinical studies and diagnosis 
of SARS-CoV-infected patients, on the characterization of 
SARS-CoV structural [4,5] and nonstructural [6,7] proteins, 
as well as on vaccines targeting the spike (S) protein [8—10]. 
Li et al. [11] discovered that angiotensin-converting enzyme 
2 (ACE2) binds efficiently to the S1 domain of the SARS- 
CoV spike protein, thus making it likely to act as the func- 
tional receptor for SARS-CoV. This is further corroborated 
by the observation that transfected NIH-3T3 cells expressing 
ACE2 can be infected with SARS pseudovirus [12]. Further- 
more, human monoclonal antibodies against S1 protein block 
the ACE2 receptor association, and enable potent neutraliza- 
tion of SARS-CoV [10]. 

Intriguingly, in stark contrast to the two known human 
coronaviruses HCoV-299E and HCoV-OC43 which usually 
cause common cold-like symptoms, SARS-CoV exerts a more 
devastating clinical effect. Host—virus interactions involve an 
elaborate interplay between molecular pathways initiated by 
the host in response to viral replication, counteracted by 
mechanisms triggered by the virus to favor its propagation. 
Insights into the host antiviral strategies as well as viral patho- 
genesis can be achieved through the elucidation of these host- 
and virus-mediated transcriptional responses, which may pro- 
vide clues to disease manifestations caused by the virus. 

Here we report the application of cDNA microarray analy- 
sis of the changes in transcriptional activities of African green 
monkey (Cercopithecus aethiops) kidney clone E6 (Vero E6) 
cells induced by SARS-CoV infection. Unlike many animal 
coronaviruses, which infect only cells of their natural host 
species, SARS-CoV can infect and replicate rapidly in cul- 
tured cells of different but related species, such as Vero 
E6 cells, within a 24-h period [13]. In addition, expression 
profiles of specific genes, including differentially expressed 
in normal and neoplastic cells or MAP-kinase activating death 
domain (DENN/MADD), Fas apoptotic inhibitory molecule 
(FAIM), inhibitor of apoptosis protein | (hIAP1) and myx- 
ovirus (influenza virus) resistance | (Mx1), were analyzed. 
We also investigated expression of the p53 tumor suppressor 
gene and its alternative splice variant in SARS-CoV-infected 
cells. These genes play important roles in apoptosis and in 
host antiviral resistance. Our study can serve as a basis for 
more detailed investigations into the exact roles of critical 
genes and pathways in the pathogenesis of this newly emerg- 
ing lethal infectious disease of humans. 


2. Materials and methods 


2.1. Cell culture and virus infection 


Vero E6 cells (ATCC number CRL-1586) were cultured 
in M199 medium (Invitrogen, Carlsbad, CA) supplemented 
with 10% fetal calf serum, 2.2 g I! sodium bicarbonate, 
10 mM HEPES, 0.1 mM non-essential amino acids, 1 mM 
sodium pyruvate at 37 °C under 5% CO, in a humidified incu- 
bator. When cells reached ~90% confluency, one batch of Vero 


E6 cells that served as the control was mock infected with 
sterile medium [14,15]. Another batch was infected with 
SARS-CoV (2003 VA2774 isolated from a SARS patient in 
Singapore) at a multiplicity of infection (MOD) of 1 [13], with 
a virus inoculum volume of 0.75 ml diluted with 1.25 ml of 
maintenance medium with only 3% fetal calf serum. After 
adsorption for | h, the inoculum was removed, and 20 ml of 
maintenance medium was added. At 12 h post-infection, cyto- 
pathic effects (CPE) were observed in infected cells. Follow- 
ing incubation at 37 °C for 12 h, both uninfected and infected 
Vero E6 cells were harvested for RNA extraction. This pro- 
tocol was designed to minimize any effects on gene expres- 
sion from “carry-over” cytokines or chemokines in the cell 
culture supernatant. Moreover, Vero cells are deficient in the 
interferon (IFN)-alpha/beta response [16]. Viral culture and 
RNA isolation were carried out in a biosafety level 3 labora- 
tory. 


2.2. Total RNA preparation 


Total cellular RNA was extracted from uninfected and 
SARS-CoV-infected Vero E6 cells using an RNeasy Mini 
Total RNA extraction kit (Qiagen, Hilden, Germany). The 
RNA was suspended in nuclease-free water, and the concen- 
tration quantitated by UV spectrophotometry at 260 nm. The 
integrity of the extracted RNA was confirmed by electrophore- 
sis ina 1% denaturing agarose gel. 


2.3. Microarray hybridization and analysis 


Total RNA (25 ug) extracted from each batch of unin- 
fected and SARS-CoV-infected Vero E6 cells was reverse tran- 
scribed and labeled with Cy3- and Cy5-dCTP, respectively, 
using a CyScribe first-strand cDNA labeling kit (Amersham 
Biosciences, Piscataway, NJ) with anchored oligo(dT) and 
random nonamers. Labeled cDNAs were purified using a 
CyScribe GFX purification kit (Amersham Biosciences) and 
resuspended in 60 ul of elution buffer. The amounts of Cy3- 
and Cy5-labeled cDNAs were measured at absorbances of 
550 and 650 nm, respectively, and calculated based on the 
following formula: (absorbance/extinction coefficient) x vol- 
ume of cDNA x dilution factor x 10'*. The extinction coef- 
ficients for Cy3 and Cy5 were 150,000 and 250,000 1 mol"! 
cm’, respectively. 

The experiments were performed using the IntelliGene II 
Human 16K cDNA microarray set (Takara Bio, Shiga, Japan), 
which comprises a total of 16,000 human genes spotted onto 
four slides, each consisting of 4000 genes. This cCDNA-based 
array consists of ~200—1000-bp fragments representing a wide 
range of known genes, hypothetical genes and expressed 
sequence tags (ESTs). 

The Cy3- and Cy5-labeled cDNA samples (40 pmol each) 
were hybridized to each microarray slide at 40 °C for 17 h 
within a Lucidea SlidePro hybridization chamber (Amer- 
sham Biosciences). The slides were washed according to the 
manufacturer’s recommendations. Each slide was scanned 
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with the GenePix 4000B image scanner, and the data were 
collected with the GenePix Pro 4.1 software (Axon Instru- 
ments, Union City, CA). The raw data were analyzed using 
Acuity software version 3.1 (Axon Instruments), as well as 
MIDAS software version 2.17, available at the TIGR website 
(http://www.tigr.org/software/tm4/) [17]. The typical flow of 
the MIDAS involved the low-intensity filter, LOWESS nor- 
malization, standard deviation regularization, and SLICE 
analysis [18]. 


2.4. cDNA synthesis and real-time SYBR Green RT-PCR 
detection 


Total RNA (4 ug) was reverse transcribed in a 40-ul reac- 
tion mix containing | x first-strand buffer, 10 mM dithiothre1- 
tol, 600 ng random hexamers, 0.5 mM deoxyribonucleoside 
triphosphates, 10 U RNase inhibitor and 200 U SuperScript 
reverse transcriptase (Invitrogen). Synthesized cDNA was ali- 
quoted and stored at —80 °C. 

The differential expression of 17 genes selected after the 
microatray analysis was further validated by real-time quan- 
titative RT-PCR using SYBR Green-based detection with a 
LightCycler system (Roche, Basel, Switzerland). In addi- 
tion, the expression profiles of DENN, FAIM, hIAP1 and 
Mx! genes were included in the real-time RT-PCR analysis. 
All the primers were designed based on the corresponding 
human mRNA sequences (Table 1). SARS-CoV RNA was 
also amplified using the primers described by Peiris et al. [2]. 
Each 10-ul reaction included 1-2 ul of cDNA, 0.5 uM of each 


Table | 


primer, 2.5—-3 mM MgCl,, and | ul of ready-to-use 10 x Hot 
Start reaction mix (Roche). The reactions were subjected to 
an initial denaturation of 95 °C for 10 min, followed by 
40 cycles each of 95 °C for 10 s, 58-68 °C for 5 s and 72 °C 
for 10 s, before being subjected to melting curve analysis. 
Amplified products were also analyzed for specificity by aga- 
rose gel electrophoresis, and were further verified by auto- 
mated cycle sequencing. The conserved HUEL gene was 
selected as the housekeeping gene transcript to normalize the 
samples, as its expression level was previously confirmed to 
be reliably constant in a number of representative tissues 
[19,20]. To ensure consistency in threshold cycle (C,) val- 
ues, duplicate reactions were performed (i.e. from each of 
the two cDNA preparations of each RNA sample), and the 
mean C, values were used for calculating the relative expres- 
sion levels. The C; values were analyzed as described previ- 
ously [21,22], and the normalized Cy values of each gene 
were subjected to Student’s t-test with two-tailed distribution 
to determine the significance at the 95% confidence level. 


2.5. Classical RT-PCR amplification of p53 splice variant 


Each cDNA sample was subjected to classical RT-PCR 
amplification of the p53 splice variant using primers P53U 
and P53D spanning exons 7—10, and the products were elec- 
trophoresed in an agarose gel [23]. The target cDNA frag- 
ment size of the alternative p53 transcript was expected to be 
133 bp larger than the wild-type p53 transcript. Amplified 
fragments were cloned into pGEM T-Easy vector (Promega, 
Madison, WI) and directly sequenced. 


Sequences of primers for the amplification of selected genes evaluated by real-time RT-PCR 


Gene Primer sequence (5’ > 3’) Target size (bp) 
Sense Antisense 
ANXA2 GTGACTTCCGCAAGCTGATG CTGAATGCACTGAACCAGGT 299 
CLK1 CCAACCATGTGATGTCTGGA AGGTCAAAGAGACGCTCATG 288 
CTSL AGGCAGGTGATGAATGGCTT ACCATTGCAGCCTTCATTGC 276 
DENN* GCTAGTGGATCTGGACAG CATCCACAGAATCCACATCG 275 
DUSPI1 AGCAGAGGCGAAGCATCATC GGATGTGAAGAGCCTCACCT 270 
EEFIG TCAGACCTTCATGAGCTGCA TACTCTCGAACCAGCGTCTG 249 
FAIM* ACTACATCAGGCAAACGAGT GCACCATACGTCCATAGCAT 294 
HAVCRI ACCATGAACCAGTAGCCACT GCAAGAAGCACCAAGACAGA 265 
hIAP1* CAGAAGACACAGACGTCTTTA CGAACTGTACCCTTGATTGTA 350 
HSPAIA AAGCAGACGCAGATCTTCAC CTCGATCTCCTCCTTGCTCA 282 
HSPCA GGCAGAGGCTGATAAGAACG TCTTCCATGCGTGATGTGTC 243 
HUEL” TCAGACGACGAAGTCCCCATGAAG TCCTTACGCAATTTTTTCTCTCTGGC 130 
IGFBP3 ACAGCATGCAGAGCAAGTAGAC CTGCTCTTTGCTGACTACTGGA 245 
IL8RA CTAATTAGCATGGCCACATCT CAGACACTGCAACACACCTGA 220 
KLF5 CCACCTGTACCAGCTACTGA CTGGAGCATCTCTGCTTGTC 283 
MADH2 GCTCTTCTGGCTCAGTCTGT TAGGGACCACACACAATGCT 305 
Mx1* GCCAGGACCAGGTATACAG GCCTGCGTCAGCCGTGC 361 
NFKBIA ACCTGGTGTCACTCCTGTTG AGCTCGTCCTCTGTGAACTC 217 
NPC2 AGCTCTGCTGCTTCAACAAC AGGTGTAGAAAGAGGCCACA 223 
RAC1 CTGACCAGCTTTTGCGGAGA AGGCATGGCAGGTGTAAGAG 289 
SARS-CoV* TACACACCTCAGCGTTG CACGAACGTGACGAAT 182 
STAM1 TGGTCAGTTCCGTTCAAGGA TAGCAGGTATCTGCCACCAA 268 
TIMP2 AAGCGGTCAGTGAGAAGGAA GGCACGATGAAGTCACAGAG 253 


* Genes included for real-time RT-PCR in addition to those selected from the microarray analysis. 
> Housekeeping gene transcript serving as the normalization control for real-time RT-PCR. 
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3. Results 


3.1. Microarray analysis of SARS-CoV-induced CPE 
of Vero E6 cells reveals novel differential human gene 
expression profiles 


By 12 h post-infection, ~50% of the cell population is 
infected, ~25% of the cells round off, with resultant produc- 
tion of an enormous number of virus particles, indicating the 
suitability of using Vero E6 as host cells for SARS-CoV infec- 
tion [13]. 

The slides were scanned for Cy3 and Cy5 signals with a 
dual-color image scanner. The majority of signals were of 
relatively low intensity, which could likely be due to the 
hybridization of primate cDNAs against a human gene-based 
microarray. However, a significant number of spots dis- 
played strong intensities, thus facilitating the analysis. The 
latter spots represent highly conserved primate cDNAs, as 
evidenced by the sequence homologies of greater than 95% 
between primate and human cDNAs of most of the genes ana- 
lyzed. Stringency was exercised by flagging individual spots 
of poor quality (e.g. extreme unevenness in signals and 
streaks), which were omitted from Acuity and MIDAS soft- 
ware analyses. Under the statistical analysis program, we ini- 
tially conducted a low-intensity filtering based on 30% of the 
standard deviation of the total background intensity follow- 
ing MIDAS software analysis. We selected this criterion as it 
allowed us to effectively disregard the low-intensity spots 
while retaining those of sufficient quality to be analyzed. The 
filtered data were normalized against the Cy3-labeled unin- 
fected Vero E6 cellular cDNAs as the reference by using 
locally weighted linear regression (LOWESS) analysis, pro- 
posed as a method to remove intensity-dependent effects com- 
monly found in log,(ratio) values [24]. We employed 
log,(ratio) rather than log,)(ratio) due to the former’s ratio- 
nality of treating all regulated genes in a similar fashion, and 
due to the generation of a continuous spectrum of values. The 
next step was to regularize the standard deviations of the nor- 
malized data to adjust for the variance of the measured 
log,(ratio) values to be similar throughout the entire array, 
thereby minimizing the variability of the results. Finally, 
through SLICE analysis, we defined an intensity-dependent 
Z-score threshold to be greater than 1.96 (i.e. equivalent to 
those more than 1.96 standard deviations from the local mean), 
which permitted us to identify 70 genes that were differen- 
tially expressed at the 95% confidence level. Of these, 29 tran- 
scripts exhibited an increase, while 41 demonstrated a 
decrease in expression level greater than fourfold, based on 
MIDAS analysis. Expression data generated from both the 
MIDAS and Acuity software were generally in agreement, 
thus supporting the reliability of the microarray results. 
Table 2 groups the transcripts that were significantly altered 
at 12 h post-infection according to their known functions, 
including mRNAs encoding four hypothetical proteins and 
eight ESTs with unknown roles. 

Genes with altered transcriptional patterns belonged to a 
wide range of functional classes. For example, those involved 


in the host translational machinery, and mRNAs encoding 
structural proteins such as cytoskeleton-associated proteins 
were notably downregulated. In contrast, mRNAs of heat 
shock proteins (HSPs) were significantly elevated, a com- 
mon phenomenon in cells experiencing external stresses or 
stimuli. However, transcripts of other functional categories 
exhibited variable responses to SARS-CoV infection. 


3.2. Real-time RT-PCR analysis confirms modified 
transcription of specific genes 


To verify the results of the microarray analysis of SARS- 
CoV-infected Vero E6 cells, real-time RT-PCR was per- 
formed for 17 interesting genes whose expression levels were 
either upregulated or downregulated. Only changes in rela- 
tive expression levels greater than twofold were considered 
significant. With the HUEL gene serving as a normalization 
control, the changes in expression and the P-values are listed 
in Table 2. Our real-time analysis of HUEL gene expression 
also revealed that the difference in Cy between infected and 
uninfected samples was negligible (less than 0.5), thereby jus- 
tifying the use of HUEL for normalization (data not shown). 
Overall, the gene expression trends by real-time RT-PCR con- 
curred with the microarray data for all genes except one, 
although the magnitude of changes differed. The only discor- 
dant gene was MADH2, which mediates TGF-f signaling to 
regulate cell growth and differentiation. Expression of this 
transcript in infected cells increased by 4.65-fold by microar- 
ray analysis, but by only 1.6-fold via real-time RT-PCR. The 
differences in relative expression by these two techniques may 
be attributed to their differences in kinetics and sensitivity. In 
general, relative expression changes by real-time RT-PCR 
were greater than those by microarray analysis, especially 
for genes that were markedly overexpressed. This may be 
explained by the better sensitivity and reliability of real-time 
RT-PCR compared to microarray, whose drawback is the satu- 
ration of fluorescent signals. 

The expression profiles in SARS-CoV-infected cells of 
additional genes absent from the microarray (DENN, FAIM, 
hIAP1 and Mx1) were also elucidated by real-time RT-PCR 
(Table 2). DENN, FAIM and hIAP1 participate in the anti- 
apoptotic mechanisms of host cells. Transcripts of DENN and 
hIAP1 were upregulated 3.75- and 3.54-fold, respectively, 
whereas FAIM was downregulated 5.46-fold. The mRNA of 
Mx1, an antiviral defense gene, was highly upregulated 59.5 1- 
fold. Real-time RT-PCR targeting a 182-bp fragment within 
the SARS-CoV RNA-dependent RNA polymerase gene pro- 
vided an indication of the intensity of virus replication and 
viral load. Twelve hours after infection, we observed a 9.42 x 
10*-fold rise in SARS-CoV transcripts, corresponding to a 
Cy value of 5.6, reflecting the very rapid SARS-CoV replica- 
tion process within Vero E6 cells. 


3.3. The p53 splice variant is expressed in 
SARS-CoV-infected cells 


A p53 gene fragment spanning exons 7—10, from both unin- 
fected and infected Vero E6 cells, was amplified, the sequenc- 
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Table 2 
Categories of human and viral genes displaying altered transcription profiles during SARS-CoV infection of Vero cells 
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GenBank accession 


Fold change in transcription 


no. Gene Acuity MIDAS Real-time RT-PCR 
Cell cycle and development 
NM_004417 Dual specificity phosphatase | (DUSP1) 19.15 11.12 42.67 (0.012) 
NM_000136 Fanconi anemia, complementation group C (FANCC) TAl 4.64 
NM_012206 Hepatitis A virus cellular receptor 1 (HAVCR1) -11.11 -11.18 —10.27 (0.030) 
NM_000598 Insulin-like growth factor binding protein 3 (IGFBP3) 6.67 8.72 —16.62 (0.001) 
NM_005901 MAD, mothers against decapentaplegic homolog 2 (Drosophila) (MADH2) 5.71 4.65 1.60 (0.009) 
NM_000582 Secreted phosphoprotein | (osteopontin, bone sialoprotein I, early T-lymphocyte -3.05 5.88 
activation |) (SPP1) 
Cellular metabolism 
NM_001912 Cathepsin L (CTSL), transcript variant | —2.54 4.36 4.96 (0.017) 
NM_000127 Exostoses (multiple) 1 (EXT1) 2.86 4.13 
NM_002410 Mannosy] (alpha-1,6-)-glycoprotein beta-1,6-N-acetyl-glucosaminyltransferase 11.18 7.28 
(MGATS) 
Cytoskeleton-associated proteins 
NM_018477 Actin-related protein 10 homolog (Saccharomyces cerevisiae) (ACTR10) -6.41 8.62 
NM_003388 Cytoplasmic linker 2 (CYLN2) 3.33 4.50 
NM_015282 Cytoplasmic linker associated protein 1 (CLASP1) 6.02 -7.17 
NM_021103 Thymosin, beta 10 (TMSB10) —2.63 4.09 
NM_021109 Thymosin, beta 4, X-linked (TMSB4X) 4.44 -5.92 
DNA repair 
AY275681 MutS homolog 3 (Escherichia coli) (MSH3) 7.63 6.87 
Heat shock and stress response 
NM_006145 DnaJ (Hsp40) homolog, subfamily B, member 1 (DNAJB1) 11.02 8.56 
NM_006644 Heat shock 105/110 kDa protein 1 (HSPH1) 15.61 9.40 
NM_005345 Heat shock 70 kDa protein 1A (HSPA1A) 14.37 7Ad 54.19 (0.002) 
NM_005348 Heat shock 90-kDa protein 1, alpha (HSPCA) 7.00 4.06 7.09 (0.014) 
mRNA processing factors 
NM_004593 Splicing factor, arginine/serine-rich 10 (transformer 2 homolog, Drosophila) 4.38 5.09 
(SFRS 10) 
Protein modulators 
NM_004071 CDC-like kinase 1 (CLK1) 13.13 8.09 58.49 (0.005) 
NM_007067 MYST histone acetyltransferase 2 (MYST2) 4,22 4.07 
NM_005742 Protein disulfide isomerase-related protein (P5) 5.95 5.80 
NM_002714 Protein phosphatase 1, regulatory subunit 10 (PPP1R10) 8.98 8.60 
NM_003255 Tissue inhibitor of metalloproteinase 2 (TIMP2) -8.00 -9.69 4.77 (0.012) 
Signal transduction 
NM_004039 Annexin A2 (ANXA2) -3.05 4.39 —4.04 (0.009) 
NM_000634 Interleukin 8 receptor, alpha (IL8RA) 4.44 4.47 26.91 (0.002) 
NM_006860 RAB, member of RAS oncogene family-like 4 (RABL4) 9.91 8.99 
NM_006908 Ras-related C3 botulinum toxin substrate | (rho family, small GTP-binding pro- 4.35 4.43 3.15 (0.020) 
tein Racl) (RAC1) 
NM_001040 Sex hormone-binding globulin (SHBG) 5.80 5.77 
NM_003473 Signal transducing adaptor molecule 1 (STAM1) 6.56 4.29 13.74 (0.023) 
Trafficking and protein degradation 
NM_014713 Lysosomal-associated protein transmembrane 4 alpha (LAPTM4A) 4.59 -5.10 
NM_006432 Niemann-Pick disease, type C2 (NPC2) -7.81 8.93 35.38 (0.013) 
NM_003900 Sequestosome | (SQSTM1) 6.20 4.90 
Transcription and transcriptional regulators 
NM_000937 Polymerase (RNA) II (DNA directed) polypeptide A, 220 kDa (POLR2A) 9.24 6.18 
NM_181552 Cut-like 1, CCAAT displacement protein (Drosophila) (CUTL1) 2.51 4.54 
L17131 High-mobility group protein (HMG-I(Y)) gene 38.46 -13.70 
NM_002158 Human T-cell leukemia virus enhancer factor (HTLF) 7.02 5.74 
NM_001730 Kruppel-like factor 5 (intestinal) (KLF5) 7.04 4.74 44,94 (0.001) 
NM_005594 Nascent-polypeptide-associated complex alpha polypeptide (NACA) -10.75 -10.75 
NM_020529 Nuclear factor of kappa light polypeptide gene enhancer in B-cells inhibitor, alpha 11.50 5.88 61.18 (0.008) 


(NFKBIA) 


(continued on next page) 
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Table 2 
(continued) 


GenBank accession Gene 


Fold change in transcription 


no. Acuity MIDAS Real-time RT-PCR 

NM_003457 Zinc finger protein 207 (ZNF207) 4.44 4.48 

Translation 

NM_013417 Isoleucine-tRNA synthetase (IARS) 8.13 -6.91 

NM_033251 Ribosomal protein L13 (RPL13) —5.05 -7.16 

NM_000979 Ribosomal protein L18 (RPL18) 4.35 -5.92 

NM_000980 Ribosomal protein L18a (RPL18A) -7.94 7.61 

NM_000990 Ribosomal protein L27a (RPL27A) 4.59 —§.25 

NM_000967 Ribosomal protein L3 (RPL3) 2.96 -§.23 

NM_033643 Ribosomal protein L36 (RPL36) 8.55 -13.75 

NM_001014 Ribosomal protein S10 (RPS10) -3.51 4.56 

NM_001018 Ribosomal protein S15 (RPS15) 4.23 5.38 

NM_022551 Ribosomal protein S18 (RPS18) 8.47 8.59 

NM_001022 Ribosomal protein S19 (RPS19) -7.81 -7.90 

NM_002952 Ribosomal protein S2 (RPS2) 3.28 4.76 

NM_001005 Ribosomal protein S3 (RPS3) 2.39 4.28 

NM_001402 Translation elongation factor | alpha | (EEF1A1) 4.17 -5.77 

NM_001404 Translation elongation factor | gamma (EEF1G) -3.60 4.13 4.13 (0.004) 

NM_017646 tRNA isopentenyltransferase 1 (TRIT1) -3.80 4.38 

Hypothetical proteins 

BC047078 Hypothetical protein LOC283507 4.13 8.49 

NM_024511 Hypothetical protein MGC4701 5.54 5.03 

BC040148 Hypothetical protein MGC49942 3.48 4.46 

D63478 KIAA0144 gene 4.74 4.80 

Expressed sequence tags 

AK123640 cDNA FLJ41646 fis, clone FEBRA2024019 5.74 441 

AA968926 Op38b10.s1 10.41 10.75 

AI379149 Tc66a01.x1 4,33 -5.40 

R60834 Yh04d12.1r1 7.72 7.50 

N94144 Za25b09.r1 7.21 TAI 

W91949 Zh47e08.r1 6.36 6.86 

AA677283 Zj61d05.s1 4.88 4.76 

AA476783 Zw94h11.s1 2.88 6.93 

Additional genes 

U44953 DENN mRNA N.A. N.A. 3.75 (0.017) 

NM_018147 Fas apoptotic inhibitory molecule (FAIM) N.A. N.A. 5.46 (0.013) 

AFO070674 Inhibitor of apoptosis protein-1 (hIAP1) N.A. N.A. 3.54 (0.010) 

NM_002462 Myxovirus (influenza virus) resistance 1, interferon-inducible protein p78 (mouse) N.A. N.A. 59.51 (0.004) 
(Mx1) 

AY268070 SARS coronavirus RNA-directed RNA polymerase N.A. N.A. 9.42 x 10* (0.009) 


The microarray data were analyzed using both Acuity and MIDAS software. For the real-time RT-PCR data, the values in brackets indicate the P-values at the 
95% confidence level. The test of significance was performed based on the differences between duplicate values as well as between data for infected versus 


uninfected samples. N.A., not applicable. 


ing of which corresponded to nucleotides 779-1131 of the 
wild-type African green monkey p53 gene (GenBank acces- 
sion no. X16384). Intriguingly, an additional larger fragment 
(Fig. 1) from SARS-CoV-infected Vero E6 cells was also 
amplified, sequencing of which revealed that it spanned nucle- 
otides 779-1095 and 1096-1131 of the p53 gene. Sequence 
analysis confirmed that this fragment contained an insertion 
of 133 nucleotides from intron 9, as previously reported in 
the MOLT-4 human acute lymphoblastic leukemia cell line 
[23]. Within this insertion is a premature termination codon 
upstream of exon 10, which is thus predicted to generate a 
truncated p53 product lacking the C-terminal amino acids 
encoded by exons 10 and 11. Compared with the genomic 
sequence of human p53 gene (GenBank accession no. 


X54156), the sequence of the 133-nucleotide insertion 
detected in infected Vero E6 cells harbored two point muta- 
tions, three base insertions and an 8-bp deletion, correspond- 
ing to that identified in MOLT-4 cells [23]. 


4. Discussion 


Fever, dry cough, dyspnea, myalgia and lymphopenia are 
common characteristic clinical features of SARS patients 
[2,25,26], although these are also symptomatic of other res- 
piratory infections. Chest radiographs of severely ill patients 
display rapid progression of unilateral peripheral air-space 
consolidation to bilateral patchy consolidation usually within 
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Fig. 1. RT-PCR analysis of p53 transcript variants. Lane M shows the DNA 
ladder markers. A 130-bp target fragment of the HUEL housekeeping gene 
transcript served as a normalization control to ensure equal amounts of star- 
ting cDNA template in uninfected (lane 1) and SARS-CoV-infected (lane 2) 
samples. The wild-type p53 gene was amplified as a 353-bp fragment in 
both uninfected (lane 3) and SARS-CoV-infected (lane 4) Vero E6 RNAs. 
However, an additional 486-bp fragment was amplified only from infected 
cells (lane 4), sequencing of which revealed its identity as the C-terminal 
alternatively spliced p53 transcript. 


a week [27]. How this novel human CoV causes such exten- 
sive damage compared with other known human CoVs that 
are common etiologic agents of relatively mild upper respi- 
ratory illness remains to be fully elucidated. Interestingly, ani- 
mal CoVs generally cause more severe diseases, and the civet 
cat is currently the prime suspect as the animal source of 
SARS-CoV that may have crossed the species barrier to 
humans. Examples of animal viruses crossing over to humans 
already exist, e.g. Nipah virus (pigs) and HSN 1 influenza virus 
(poultry). 

In order to understand the pathophysiology of SARS-CoV 
infection and host cellular responses, we employed cDNA 
microarray, RT-PCR and cDNA sequencing analyses to inves- 
tigate the gene expression changes during the infection pro- 
cess. Similar studies have previously been conducted on infec- 
tions with dengue virus [28], human polyomavirus JC [29] 
and herpes simplex virus [14]. Vero E6 served as the host 
cells for infection, as this cell line is the most permissive for 
SARS-CoV replication [1-3]. So far, human cell lines are 
not suitable, with the recent exception of LoVo cells, which 
are susceptible to SARS-CoV infection [30]. However, per- 
sistent infection occurs in LoVo cells with no visible CPE, in 
contrast to the lytic infection seen in Vero E6 cells which 
mimics the apoptosis observed in the respiratory epithelial 
cells of severely ill SARS patients [3 1-33]. The infection pro- 
cess is very rapid, such that virus particles are internalized 
into vacuoles within 10 min post-infection [34], extracellular 
virus particles exist in ~30% of the cells by 6 h post- 
infection, and numerous large vacuoles containing mature 
virus are present in the cytoplasm of infected cells from 12 to 
21h post-infection [13], as evident from the almost 10°-fold 
increase in the level of SARS-CoV transcripts by real-time 
RT-PCR analysis. Expression of the novel SARS-CoV pro- 
tein 3a is detectable in infected Vero E6 cells at 8—12 h post- 
infection, and also in the pneumocytes of a SARS patient’s 


lung [35]. Furthermore, the incubation period of SARS is vari- 
ously reported to be 6.4 days [36] or 5 days [37], with a range 
of 3-8 days [38]. Therefore, 12 h after infection of Vero cells 
represents an appropriate time-point that allows for a larger 
population of cells to be infected so that significant changes 
representative of the pathophysiological process of SARS- 
CoV infection can be elucidated. It would be ideal to conduct 
time-course experiments commencing from earlier time- 
points in the virus replication cycle. 

Vero E6 cells are derived from African green monkey, and 
the limitations of cross hybridization may affect our microar- 
ray results. It is estimated that 95% of nucleotides are exactly 
shared between human and chimpanzee DNA [39]. Thus, the 
closest relative of humans is the chimpanzee, followed by the 
gorilla, orangutan and Old World monkeys (macaque and Afri- 
can green monkey), in increasing order of evolutionary diver- 
gence [40,41]. In view of the relatively small differences 
between the genomes of humans and primates [42], with most 
disparities arising from insertions, deletions and rearrange- 
ments [43], Vero cells represent a viable alternative model 
for investigating cellular gene responses to SARS-CoV in lieu 
of a human cell line. This was further substantiated by 
sequencing of the genes selected for real-time analyses, which 
revealed sequence differences of only 1—5% (data not shown). 

From the microarray analysis, we identified 70 transcripts 
with altered expression that were classified according to their 
functional roles. The most distinguishable genes to be down- 
regulated by the infection were those involved in the host 
translational machinery, including 40S ribosomal proteins 
(RPS2, RPS3, RPS10, RPS15, RPS18, RPS19); 60S riboso- 
mal proteins (RPL3, RPL13, RPL18, RPL18A, RPL27A, 
RPL36); isoleucine tRNA synthetase (ARS) and modifying 
enzyme (TRIT1); and eukaryotic translation elongation fac- 
tors (EEF1A1 and EEF1G). Such host translational shutoff is 
frequently observed, as evident in infections with herpes sim- 
plex virus [44], influenza virus type A [45] and poliovirus 
[46]. However, unlike poliovirus infection, which often results 
in selective translation of viral uncapped mRNAs, CoVs give 
rise to mRNAs that are structurally similar to their eukary- 
otic hosts. This allows CoVs to parasitize upon the host 
machinery to translate the viral mRNA. Another study also 
indicates that murine CoV at high MOT inhibits host protein 
synthesis during the very early stage of infection, and sug- 
gests that the increased number of viral mRNAs produced 
during the later phase of infection compete with cellular 
mRNAs for cellular ribosomes [47]. 

SARS-CoV-infected cells also exhibited diminished 
expression of genes related to the maintenance of cytoskel- 
etal structure. CYLN2 is associated with Williams syn- 
drome, a multisystem developmental disorder [48]. It medi- 
ates the interaction between specific membranous organelles 
and microtubules, and may be an anti-catastrophic factor [49]. 
CLASP! binds CYLN proteins and microtubules, colocaliz- 
ing at the distal ends to provide stabilizing effects on the 
microtubules [50]. The microtubule network plays an impor- 
tant role in viral replication and viral protein trafficking [51]. 
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Thus, disrupting the network may initially result in less virus 
production until a point at which the cells are overwhelmed 
by virus multiplication and then lyse to release virus par- 
ticles. By exposing peritoneal macrophages to cytochalasin 
B, murine CoV infection changes from an acute cytopathol- 
ogy to a persistent type [52]. In SARS-CoV infection, these 
microtubule networks may become disrupted, making it con- 
ducive for persistent viral infection and subsequent release. 
Besides microtubules, the actin assembly may also be greatly 
affected by downregulation of TMSB4X and TMSB 10, which 
are actin-sequestering proteins. Indeed, actin mRNAs are 
reduced in murine CoV infection [47]. Furthermore, reduc- 
tion of TMSB4X transcripts may disrupt cellular functions, 
since TMSB4X may have a unique integrative function that 
links the actin cytoskeleton to important immune and cell 
growth-signaling cascades [53]. 

The upregulation of several HSPs was not unexpected 
given that elevation of HSP expression represents a crucial 
response towards external stimuli or stresses such as infec- 
tions. Noteworthy is the strong induction of HSPA1A which 
is widely involved in translocation of membrane proteins, and 
acts as a scavenger of degraded peptides for antigen presen- 
tation. During a viral infection, a rise in the level of secreted 
cytokines such as interleukin 2 (IL-2) is often encountered. 
Besides inducing the proliferation of T lymphocytes, IL-2 also 
induces the expression of HSPA1A and HSPCA [54]. In cer- 
tain inflammatory lung diseases with marked accumulation 
of eosinophils in the bronchoalveolar lavage, free oxygen radi- 
cals are produced which induce alveolar macrophages to syn- 
thesize antioxidants such as HSPs [55]. Such a phenomenon 
may operate in the atypical pneumonia characteristic of SARS 
patients who show significantly elevated plasma levels of IFN- 
gamma, IL-1, IL-6 and IL-12, IL-8, monocyte chemoattrac- 
tant protein-1, and IFN-gamma-inducible protein-10 [56]. 

In our study, elevated expression of MGATS may contrib- 
ute to the initiation of GlcNAc betal,6 branching on 
N-glycans, thereby increasing N-acetyllactosamine, the ligand 
for galectins. In turn, galectins modulate T-cell proliferation 
as well as enhancing apoptosis at the site of antigen presen- 
tation [57]. Cathepsin L (CTSL) is a lysosomal cysteine pro- 
teinase that plays a major role in intracellular protein catabo- 
lism. Suppression of CTSL expression in A549 lung cancer 
cells leads to growth inhibition, which is compatible with 
CTSL downregulation during SARS-CoV infection. How- 
ever, this growth inhibition is partially compensated by 
upregulation of IL-8 production [58]. In addition, CTSL 
secreted from human fibroblasts in response to external stimuli 
plays an important role in the processing of IL-8 to mature 
form in inflammatory sites [59]. 

Altered transcription of signal transduction molecules was 
also observed in SARS-CoV-infected cells. ANXA2, a fibrin- 
olytic receptor, binds plasminogen and tissue plasminogen 
activator (t-PA) independently at the cell surface of mono- 
cytes and macrophages, thereby enhancing the catalytic effi- 
ciency of plasmin production [60], and limiting pulmonary 
fibrosis. However, the downregulation of ANXA2 in SARS- 


CoV infection is expected to decrease fibrinolysis. IL-8 is a 
member of the family of pro-inflammatory cytokines with 
roles in chemotaxis and activation of monocytes, selective 
chemotaxis of memory T cells, and induction of neutrophil 
infiltration in vivo. In dengue virus infection, IL-8 transcrip- 
tion is controlled by NF«B [61], which activates many immu- 
noregulatory genes in response to pro-inflammatory stimuli. 
In our analysis, the upregulation of IL8RA (a low-affinity 
receptor for IL-8) implies more active IL-8 signal transduc- 
tion at the site of inflammation, where the relative concentra- 
tion of IL-8 is high [62]. 

The STAM! protein contains an SH3 domain and an immu- 
noreceptor tyrosine-based activation motif (ITAM). 
STAM1 associates with JAK3 and JAK2 kinases via its ITAM 
region, is phosphorylated by the JAK kinases upon cytokine 
stimulation, thus acting as an adaptor molecule involved in 
the downstream signaling of cytokine receptors [63]. 
STAM1 upregulation in our study suggests its enhanced role 
in the downstream signaling of cytokine receptors during 
SARS-CoV infection. In our model, several transcriptional 
modulators with altered expression may be crucial to the 
pathogenic outcome of infection. 

The TIMP gene family encodes natural inhibitors of matrix 
metalloproteinases (a group of peptidases involved in degra- 
dation of the extracellular matrix) that regulate tumor growth, 
progression, and angiogenesis. TIMP2 activates Ras through 
a PKA-mediated pathway, leading to the formation of the 
Ras/PI3K complex, with roles in proliferation, differentia- 
tion, membrane ruffling, and prevention of apoptosis [64]. 
Reduced TIMP2 expression in SARS-CoV infection is thus 
predicted to promote apoptosis of infected cells. 

Certain transcripts involved in cell cycle and development 
(e.g. DUSP1 and FANCC) were upregulated. DUSP1 plays 
an important role as a protein phosphatase in human cellular 
response to environmental stress, and as a negative regulator 
of cellular proliferation by inactivating mitogen-activated pro- 
tein kinases (MAPK) [65]. IGFBP3 is essential for TNFa- 
induced apoptosis [66], but downregulation of its transcript 
level in SARS-CoV-infected cells may imply a mechanism 
to avert virus-induced apoptosis, thus favoring viral replica- 
tion. 

We also analyzed the transcriptional profiles of several 
apoptosis-related genes. SARS-CoV-infected Vero E6 cells 
exhibited upregulation of the combined transcripts of the anti- 
apoptotic DENN isoform and the pro-apoptotic IG20 iso- 
form that both interact with tumor necrosis factor receptor 
1 and participate in the MAPK pathway [67,68]. Increased 
transcript level of the hIAP1 apoptotic inhibitor [69] during 
SARS-CoV infection mirrored that observed in dengue virus 
infection of human umbilical vein endothelial cells [28]. 
Although infected Vero cells eventually lyse at the late phase 
of infection, hIAP1 may attenuate apoptosis at the early infec- 
tion phase, thereby allowing the virus to replicate to high titer. 
In contrast, another anti-apoptotic gene, FAIM, was down- 
regulated in SARS-CoV infection. The susceptibility of pri- 
mary B cells to Fas-mediated apoptosis is modulated by FAIM 
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[70] and NF«B [71,72]. Thus, the lower expression of FAIM 
and the higher expression of NFKBIA which inhibits NF«kB 
are likely to favor Fas-mediated apoptosis. In SARS-CoV 
infection, the elevated expression of NFKBIA, an inhibitor 
of the NF«B complex, may attenuate the effect of NF«B [73]. 
A critical regulator of apoptosis, NF«B is most commonly 
involved in the suppression of programmed cell death via the 
transactivation of anti-apoptotic gene expression [74]. NF«kB 
also plays a key role in controlling cell proliferation, as evi- 
denced by the deregulated expression of its constitutively 
active form in certain cancers such as Hodgkin’s disease 
[75,76]. The interplay between, and altered expression of, 
apoptosis-related genes in SARS-CoV infection may help to 
achieve viral replication before cell death. 

Mx proteins comprise a group of antiviral GTPases that 
play an important role in IFN-induced antiviral defenses [77]. 
They also act to sequester viral nucleocapsids and limit their 
accessibility for viral replication [78]. The marked elevation 
of Mx! expression (~60-fold increase) at 12 h post-infection 
suggests a highly notable antiviral response to SARS-CoV 
pathogenesis. 

In SARS-CoV-infected Vero cells, an additional larger tran- 
script of the p53 tumor suppressor gene was detected, arising 
from the use of an alternative 5’ donor splicing site, with the 
resultant loss and modification of the C-terminus [23]. Vari- 
ous alternatively spliced p53 isoforms are expressed at rela- 
tively low levels and tend to be restricted to particular cell 
types and/or physiologic conditions [79]. The p53 C-termi- 
nus plays a pivotal role in regulating the activity of the wild- 
type molecule. Compared to the regularly spliced form of 
p53, the C-terminally altered p53 protein inhibits both p53- 
dependent apoptosis and transactivation [80]; binds more effi- 
ciently to DNA in a sequence-specific manner; is more effi- 
cient in concentration-dependent transcriptional repression of 
the promoter of the p21 cyclin-dependent kinase inhibitor 
gene; associates with and interferes more efficiently with bind- 
ing of TATA-binding protein to a TATA-containing DNA 
sequence [81]. The presence of this p53 isoform with a trun- 
cated C-terminus is therefore expected to modify the activity 
of the wild-type p53, and exert a role in the pathogenesis of 
SARS-CoV infection. 

Our study offers an overview of the cascade of changes in 
host cellular expression culminating from infection with 
SARS-CoV. The counterbalancing of several anti- 
inflammatory and pro-inflammatory pathways together with 
the variable expression of apoptosis-related genes was strik- 
ing and instructive [82]. These molecular events underlie the 
mechanism that ensures cell survival during the early phase 
of infection to allow rapid multiplication of progeny virus 
before the CPE occur. Being deficient in IFN response, Vero 
cells may not necessarily represent physiologically normal 
cells. Nonetheless, alterations in transcription of genes that 
are conserved between African green monkeys and humans 
provide a better understanding of the biology of SARS-CoV 
infection. Such molecular insights into the pathophysiologi- 
cal mechanisms may lead to novel and viable strategies for 
intervention in the infection process [83]. 
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